Core Precession and Global Modes in Granular Bulk Flow 
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A transition from local to global shear zones is reported for granular flows in a modified Couette 
cell. The experimental geometry is a slowly rotating drum which has a stationary disc of radius 
Rs fixed at its bottom. Granular material, which fills this cell up to height H, forms a wide 
shear zone which emanates from the discontinuity at the stationary discs edge. For shallow layers 
{H/Ra < 0.55), the shear zone reaches the free surface, with the core of the material resting on the 
disc and remaining stationary. In contrast, for deep layers (H/Rs ^ 0.55), the shear zones meet 
below the surface and the core starts to precess. A change in the symmetry of the surface velocities 
reveals that this behavior is associated with a transition from a local to a global shear mode. 

PACS numbers: 45.70.Mg, 45.70.-n, 83.50.Ax, 83.85.Cg 



Slowly sheared granular matter usually organizes into 
rigid regions separated by narrow shear bands where the 
material yields and flows Iilllllllll3,Elli- Such 
grain flows appear to prohibit successful continuum mod- 
eling - not only because of the steep gradients in velocity, 
but also because subtle microscopic characteristics of the 
granulate can alter the flow in a qualitative manner [^. 

The formation of narrow shear bands can be avoided by 
driving the granulate from a discontinuity in the bottom 
support of the grain layer (Fig.^), effectively pinning a 
wide shear zone away from the sidewalls . The result- 
ing grain flows are smooth and robust, with both velocity 
profiles and shear zone location exhibiting simple, grain 
independent properties - these flows should be amenable 
to a continuum description. At present, a simple model 
aimed at describing an infinitely narrow shear band cap- 
tures the location of the shear zone |0| , but there are no 
models which describe the finite width or velocity profile 
of the shear zones found in experiment. In this Letter 
we uncover a novel transition of the fiow in this system 
which occurs for deep grain layers. 

The flows of shallow layers of granular materials in our 
setup are characterized by wide shear zones. The surface 
velocities are azimuthal and proportional to the driving 
rate fi, and the ratio of the average angular surface ve- 
locity and fl, denoted uj{r), is well fitted by 0: 

u}{r) ~ nerf ( —) , (1) 
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where r is the radial coordinate, Rc and W parameterize 
the location and width of the shear zones, and nerf(a;) 
denotes the normalized errorfunction 1/2-1-1/2 errorf(a;). 
W is independent of disc radius Rg and grows with grain 
diameter d and layer height H as 9] 

W/d - {H/df/^ H^/^d^/^ 

while Rc is independent of d 

1 - Rc/Rs = {H/Rs 



(2) 
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When H/Rs is small, lo tends to zero for small r 
core material remains stationary in shallow layers. 



(3) 
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FIG. 1: (a) Schematic side-view of our split-bottomed shear 
cell, showing a stationary bottom disk of radius Rs (grey), the 
rotating bottom ring (striped) and rotating outer cylinder of 
radius Ro — 105 mm. (b-f ) Series of snapshots of top- views of 
our setup (for Rs — 95 mm and -ff = 60 mm) , where colored 
particles sprinkled on the surface illustrate the core precession 
for t = s (b), t = 10 s (c), t = 10^ s (d), t = W^ s (e) and 
t = 10^ s (f). 



inward shift of Rc with H/Rs (Eq. implies, however, 
that this behavior should break down eventually. 

Figs, n^-f illustrate novel flow patterns, characteristic 
for deep layers. The most striking feature is that the 
core starts to precess with a nonzero rate ujp — Lu{r -^0). 
Precession is not simply the consequence of the overlap 
of the two opposing shear zones given by Eq. (1) at r = 
±i?c) since before being eroded by shear, the inner core 
rotates as a solid blob for an appreciable time (Fig. l(b- 
f))- 

In this Letter we will characterize the transition to pre- 
cession, and specifically address the following questions: 
How does LOp grow with layer depth? What are the ve- 
locity profiles for such deep layers? What happens to the 
scaling relations Eq. 2-3? What underlies this transition? 

Setup - Our setup is a modified version of the disk 
geometry described in earlier work , and consists of a 
stationary bottom disk of radius Rs, a rotating bottom 
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FIG. 2: (a) Surface velocity profiles Lu{r) for Rs = 95 mm 
and increasing layer depth. Thick curves: H — 10, 20, . . . , 80 
mm; Thin curves H = 15, 25, ... 75 mm; Dashed curves H = 
56, 56, ... 69 mm. Precession gradually sets in for // > 60 
mm. The thick line indicates the nonzero plateaux value Up 
for H — 65 mm. (b) Precession rate uip as function of H for 
Rs = 45 mm (diamonds) , -Rs = 65 mm (x) and Rs = 95 mm 
(circles), (c) Precession rate curves collapse when plotted as 
function of H/Rs. The inset shows that grows smoothly 
with H/Rs. 



ring and outer cylinder of radius i?o = 105 mm (Fig. la). 
The disc radius Rs can be varied from 45 mm to 95 mm. 
The cell is filled to a height H with a polydisperse mix- 
ture of spherical glass beads with diameters ranging from 
0.6 to 0.8 mm; a layer of grains is glued to the side walls 
and bottom rings to obtain rough boundaries. Slow flows 
are rate-independent, i.e., surface velocities are propor- 
tional to the rotation rate Q of outer cylinder and its 
co-moving ring. We fix at 0.15 rad/s, and recover the 
surface velocity by a variant of particle image velocimetry 
(for more details see 0). 

Basic phenomenology - Figs. illustrate the main 
features of the surface velocity profiles uj{r) as function 
of increasing layer depth H. For shallow layers {H/Rg ^ 
0.5) the data is well fitted by Eq.^ The core precession, 
which sets in for H/Rg < 0.65 constitutes the clearest 
deviation from this simple form. For H/Rs ^ 0.8 cUp 
tends to one, and the whole surface rotates then rigidly. 
Apparently the shear zone remains confined below the 
surface in this case. 

In Fig. 1213-c we show that the data for usp collapses 
when plotted as function of H/Rg. The data suggests 
a much faster than exponential growth in the initial 
stages (i.e. when ujp <C 1) as illustrated in the inset 
of Fig. 13;). Before precession becomes appreciable, the 
left-right symmetry of uj{r) is broken: uj{Ar — Rc) ^ 
1 — llj{Rc — Ar). As illustrated in Fig.Ol the right (large 
r) tail of uj{r) is then significantly steeper than its left 
(small r) tail. 




FIG. 3: Illustration of broken left-right symmetry of lu for 
Rs = 95 mm and H — 58 mm, where A denotes the rescaled 
spatial coordinate ((r — Rc)/Rs)- Shown are (jj{X) (grey), its 
symmetric counterpart 1 — a;(— A) (dashed) and the best fit of 
Lj(A) to Eq. Q (dots). 



Hence, three regimes can be distinguished. Shallow 
layers occur for H/Rg ^ 0.45 where LUp is zero, Lu{r) is 
symmetric and Eq.^fits the data well. Deep layers occur 
for H > 0.65, where the reflection symmetry of zones is 
strongly broken and precession sets in. In the crossover 
regime, 0.45 ^0.65, the symmetry of uj{r) is broken 
but LUp remains small. 

Transition from local to global modes - Both symme- 
try breaking and precession reflect the same underly- 
ing change in the qualitative nature of the shear modes. 
Since uj{r) in general can be written as nerf [x(r)], we 
have calculated x(^) '■— nerf~^ ['^('')] from our experi- 
mentally obtained velocity profiles. The resulting xi''') 
are shown in Fig.^. For shallow layers x(r) is a linear 
function, consistent with the fact that a;(r) obeys Eq. 
For deep layers, x(^) becomes a nonlinear function and 
over the whole range of H, x(r) can be fitted well by an 
empirical third order polynomial llll | (see Fig.^): 



X(r) = flo + air + a^r^ 



(4) 



Both precession and symmetry breaking are well captured 
by this fit. As shown in Fig. ^-c, the three aforemen- 
tioned regimes have now a simple interpretation in terms 
of the linear (ai) and cubic (03) coefficient. For shallow 
layers the cubic term vanishes while for deep layers the 
linear term is absent; both are nonzero in the crossover 
regime. 

Note that the fit parameters ai and 03 have dimensions 
of inverse length and inverse length cubed, and are there- 
fore expected to scale as W and W^. Rewriting Eq. |21as 
M^/(i?'/^)di/3 ~ (F/i?,)2/3, it follows that aii??/^di/3 
and a^R^d should collapse wh en p lotted as function of 
H/Rs - indeed this is the case [13 (Fig. 4c). 

The functional form for uj for deep layers is thus 

uj{r) « nerf(ao + asr^) , (5) 

while for shallow layers Eq. Q can be rewritten as 

uj{r) w ncrf(ao + air) . (6) 

These two limits have different symmetries. For shal- 
low layers, two flow profiles obtained for two different 
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H(mm) 90 0.2 H/R^ (mm) 0.8 



FIG. 4: (a) Profiles of x('') for = 95 n™i S'Hd increasing 
layer depth (dots), compared to cubic fits given by Eq. @ 
(curves). Similar to Fig. 2a, H = 10, 15, 20, ... , 55, 56, . . . , 70 
mm. (b) Fit parameters ai for Rs = 65 mm (diamonds) and 95 
mm (+) and as for Rs — 65 mm (triangles) and 95 mm (stars) 
as function of H. (c) Non-dimensionalized fit parameters as 
function of H/Rs (see text). 



slip radii (Rgi and Rs2) are simply related by transla- 
tions of the radial coordinate via w(r — Rsi)\rsi,h = 
w(r — Rs2)\rs2,h ~ the shear zone is a local phenomenon, 
insensitive to the location of the center of the shear cell. 
For deep layers this symmetry is absent, and the shear 
mode is global. The precession and symmetry breaking 
of U!{r) both reflect the crossover from a local to a global 
shear mode. 

3D flow structure - In Fig. 5, the crucial difference 
between the 3D flow structure of shallow and deep layers 
is illustrated. The flow in the bulk has been probed by 
putting patterns of lines of colored tracer particles at 
given height in the bulk, Hh, adding more material so as 
to bury the line-pattern, rotating the system for a short 
period ('-~^ 8 s), and recovering the deformed line-pattern 
by carefully removing the upper layers of grains 9J . The 
striking difference between these two cases where w 
and a Up ~ 0.8 respectively, is that in the former case 
the shear zones reach the free surface while for the latter 
case the shear zones meet in the bulk. This scenario is 
consistent with the finding of recent MRI and numerical 
studies of the flow in this geometry 0, j and gives a 
straightforward interpretation of local and a global shear 
modes. 

Comparison to theory - Recently, Unger et al. pro- 
posed a simple theory for the location of the shear zone 
in our system. They model the shear zone as an infinitely 
thin sheet along which the granulate slides, with the nor- 
mal stress on the sheet given by hydrostatic pressure and 
the shear stress by sliding friction. By numerically ad- 
justing the sheet so as to minimize the total torque, the 




FIG. 5: 3D flow-profiles for Ra = 95 mm. Both stacks show 
five slices at heights Hb = 11,21,.. .,51 mm inside the ma- 
terial; the left stack is for layer depth H = 57 mm, the right 
stack for H = 71 mm. The bottom figures show sketches of 
the shear zones within the bulk (black) - the dashed curves 
indicate the heights Hi, where the patterns where created. 

theory predicts the shear zones location at the free sur- 
face, Rm{H). This function is different from the simple 
scaling law Eq. (jSJ, and the theory predicts that for deep 
layers the shear sheet no longer reaches the free surface 
but closes in the bulk; the two cases are separated by 
a first order (hy steretic) transition that occurs around 
ff/i?, «0.7 lig. While we have not found any evidence 
for hysteresis, the theory of Unger et al. does a good job 
in capturing the qualitative behavior. 

To test the theory quantitatively, in Fig. 6 our data 
for the shear zones location is compared to Rm{H) and 
Eq. (PJ. For the general case, we have to define the center 
of general shear zones. There is no unique choice, and 
we have tested the following three expressions: 

i?i where oj{r) = 0.5 , (7) 
i?2 := where drOj{r) is maximal , (8) 




20 H(mm) 60 80 



FIG. 6: Shear zone positions versus layer depth for Rs — 95 
mm, where circles, pluses and diamonds correspond to _Ri, R2 
and -R3 respectively, the solid curve is the scaling form given 
by Eq. and the dashed curve is the result by Unger et al. 
flO|| . The inset shows the same data, now in rescaled form, 
by defining C - (1 - R^/ Rs) / {H / Rsf'^ . 
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FIG. 7: (a) Shear zones width W' as function of H for _Ra = 
45, 65 and 95 mm. For small heights the shearzone width is 
independent of R^, while deviations start to set in for larger 
heights, (c) The nondimensionalized width W' FdJ^ d}^^^ fol- 
lows a universal curve as function of H/Rs, which for small 
heights is approximated well by a powerlaw with exponent 
2/3, as indicated by the straight line. 



i?3 



where dr{rLj{r)) is maximal 



(9) 



Note that for shear zones of the form Eq. Ri and 
R2 coincide with R^, while i?3 corresponds to where the 
dimensional strain rate is maximal. As shown in Fig. El 
for shallow layers, both i?i and R2 are better described 
by Eq. |j3Jl than by i?™, while R^ appears to be better 
described by For deep layers, the situation is more 
complicated, with neither model describing any of the 
three measured curves in detail. 

Open Questions - There are three questions that we 
think deserve particular attention. The first concerns the 
functional form of While it is now widely accepted 

that in many shear flows, the velocity profile across a 
shear band is smooth, it is not clear what mechanism se- 
lects this velocity profile, or indeed its "tail" 0. In some 
systems, these tails appear exponential for our shal- 
low shear zones these tails are Gaussian l9|, and our re- 
sults here suggest that for deep layers, the tail stretching 
out into the material has a more complex form (Eq. 
Similarly, we do not know what determines the functional 
form or tails of ujp{H/Rs) (Fig. 2b-c). Theoretical work 
is very welcome here. 

The second question concerns the finite width of the 
shear zones. Let us define the width W, for general 
velocity profiles, as the interval where u grows from 
LUp + 0.1 X (1 - ujp) to 0.9 X (1 - LOp) W does not 

depend on Rg for shallow layers, but Rs sets the scale 
where W saturates; for deep layers, both the microscale 
(d) and macroscale (Rs) play a role (Fig. 7a). The ra- 
tio Rs/d between these scales can be conceived as the 
non-dimensionalized radius of curvature. Eq. (|21), when 
rewritten as W'/d = {H/ Rsf/^{Rs/df/^ , suggests to 
plot the rescaled width, W /{R^J^d^/^), as function of 
H/ Rs - as shown in (Fig. 7b) , this leads to a good data- 
collapse. This result is consistent with the scaling with 
exponent 2/3 for shallow layers (Eq. (jSJ). It is not known 
what sets the value of this non-trivial exponent. 

Finally, it is an open question whether the transition 
to precession should be conceived as a smooth crossover 



or as a sharply defined transition. The smooth growth 
of u)p with H/Rs suggests a crossover (Fig. 2b-c), while 
the (critical) vanishing of ai with (H/Rg) (Fig. 4b-c) 
suggests a sharp transition. 

Outlook - A wide range of granular flows can be 
achieved in split bottom geometries. Detailed studies of 
the 3D structure of these flows ll-dl , together with studies 
of the role of packing density anisotropies 0| and 
contact forces are uncovering the richness of slowly 
sheared granulates, and provide crucial ingredients for 
theories of flows of dense granular media. The robust be- 
havior of grain flows in split bottom geometries hopefully 
will provide important testing ground for such theories 

Mm 
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